General contraction of Gaussian basis sets. Part 2: Atomic natural orbitals and the calculation of atomic and molecular properties

A recently proposed scheme for using natural orbitals from atomic configuration interaction (CI) wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outmost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital (ANO) sets

A definitive heat of vaporization of silicon through benchmark ab initio calculations on SiF_4

In order to resolve a significant uncertainty in the heat of vaporization of silicon -- a fundamental parameter in gas-phase thermochemistry -- $\Delta H^\circ_{f,0}$[Si(g)] has been determined from a thermochemical cycle involving the precisely known experimental heats of formation of SiF_4(g) and F(g) and a benchmark calculation of the total atomization energy (TAE_0) of SiF_4 using coupled-cluster methods. Basis sets up to $[8s7p6d4f2g1h]$ on Si and $[7s6p5d4f3g2h]$ on F have been employed, and extrapolations for residual basis set incompleteness applied. The contributions of inner-shell correlation (-0.08 kcal/mol), scalar relativistic effects (-1.88 kcal/mol), atomic spin-orbit splitting (-1.97 kcal/mol), and anharmonicity in the zero-point energy (+0.04 kcal/mol) have all been explicitly accounted for. Our benchmark TAE_0=565.89 \pm 0.22 kcal/mol leads to $\Delta H^\circ_{f,0}$[Si(g)]=107.15 \pm 0.38 kcal/mol ($\Delta H^\circ_{f,298}$[Si(g)]=108.19 \pm 0.38 kcal/mol): between the JANAF/CODATA value of 106.5 \pm 1.9 kcal/mol and the revised value proposed by Grev and Schaefer [J. Chem. Phys. 97, 8389 (1992}], 108.1 \pm 0.5 kcal/mol. The revision will be relevant for future computational studies on heats of formation of silicon compounds.Comment: J. Phys. Chem. A, submitted Feb 1, 199

The determination of accurate dipole polarizabilities alpha and gamma for the noble gases

The static dipole polarizabilities alpha and gamma for the noble gases helium through xenon were determined using large flexible one-particle basis sets in conjunction with high-level treatments of electron correlation. The electron correlation methods include single and double excitation coupled-cluster theory (CCSD), an extension of CCSD that includes a perturbational estimate of connected triple excitations, CCSD(T), and second order perturbation theory (MP2). The computed alpha and gamma values are estimated to be accurate to within a few percent. Agreement with experimental data for the static hyperpolarizability gamma is good for neon and xenon, but for argon and krypton the differences are larger than the combined theoretical and experimental uncertainties. Based on our calculations, we suggest that the experimental value of gamma for argon is too low; adjusting this value would bring the experimental value of gamma for krypton into better agreement with our computed result. The MP2 values for the polarizabilities of neon, argon, krypton and zenon are in reasonabe agreement with the CCSD and CCSD(T) values, suggesting that this less expensive method may be useful in studies of polarizabilities for larger systems

A simple interpretation of the Fe2(-) photoelectron spectrum

The photoelectron spectrum of Fe2(-) can be simply interpreted in terms of electron detachment from the mildly bonding 4s sigma sup * sub u orbital of a (4s sigma sub s)(sup 2) (4s sigma sup * sub u)(sup 2) (3d)(sup 13) anion. This interpretation implies a (4s sigma sub g)(sup 2) (4s sigma sup * sub u)(sup1) (3d)(sup13) configuration for the ground state of Fe2, correlating with one ground state (4s2 3d6) and one excited state (4s1 3d7) Fe atom. A comparison of the bond length and vibrational frequency of Fe2 to values for transition metal dimers containing single 4s-4s bonds is suggestive of 3d-3d bonding in this molecule. The results of preliminary full-valence configuration interaction calculations provide stron support for the proposed Fe2 and Fe2(-) configurations

Nse2, a component of the Smc5-6 complex, is a SUMO ligase required for the response to DNA damage

The Schizosaccharomyces pombe SMC proteins Rad18 (Smc6) and Spr18 (Smc5) exist in a high-M(r) complex which also contains the non-SMC proteins Nse1, Nse2, Nse3, and Rad62. The Smc5-6 complex, which is essential for viability, is required for several aspects of DNA metabolism, including recombinational repair and maintenance of the DNA damage checkpoint. We have characterized Nse2 and show here that it is a SUMO ligase. Smc6 (Rad18) and Nse3, but not Smc5 (Spr18) or Nse1, are sumoylated in vitro in an Nse2-dependent manner, and Nse2 is itself autosumoylated, predominantly on the C-terminal part of the protein. Mutations of C195 and H197 in the Nse2 RING-finger-like motif abolish Nse2-dependent sumoylation. nse2.SA mutant cells, in which nse2.C195S-H197A is integrated as the sole copy of nse2, are viable, whereas the deletion of nse2 is lethal. Smc6 (Rad18) is sumoylated in vivo: the sumoylation level is increased upon exposure to DNA damage and is drastically reduced in the nse2.SA strain. Since nse2.SA cells are sensitive to DNA-damaging agents and to exposure to hydroxyurea, this implicates the Nse2-dependent sumoylation activity in DNA damage responses but not in the essential function of the Smc5-6 complex

On the bond distance in methane

The equilibrium bond distance in methane was optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valence correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers

Connected triple excitations in coupled-cluster calculations of hyperpolarizabilities: Neon

We have calculated the second hyperpolarizability gamma of neon using the CCSD(T) method. The accuracy of the CCSD(T) approach has been established by explicit comparison with the single, double and triple excitation coupled-cluster (CCSDT) method using extended basis sets that are known to be adequate for the description of gamma. Our best estimate for gamma(sub 0) of 110 +/- 3 a.u. is in good agreement with other recent theoretical values and with Shelton's recent experimental estimate of 108 +/- 2 a.u. Comparison of the MP2 and CCSD(T) hyperpolarizability values indicates that MP2 gives a very good description of the electron correlation contribution to gamma(sub 0). We have combined MP2 frequency-dependent corrections with the CCSD(T) gamma(sub 0) to yield values of gamma(-2 omega;omega,omega,0) and gamma(exp K)(-omega;omega,0,0)

The dissociation energy of N2

The requirements for very accurate ab initio quantum chemical prediction of dissociation energies are examined using a detailed investigation of the nitrogen molecule. Although agreement with experiment to within 1 kcal/mol is not achieved even with the most elaborate multireference CI (configuration interaction) wave functions and largest basis sets currently feasible, it is possible to obtain agreement to within about 2 kcal/mol, or 1 percent of the dissociation energy. At this level it is necessary to account for core-valence correlation effects and to include up to h-type functions in the basis. The effect of i-type functions, the use of different reference configuration spaces, and basis set superposition error were also investigated. After discussing these results, the remaining sources of error in our best calculations are examined

Wetlands and Coal Surface Mining: A Management Handbook

As the third phase of a three-year project, this report outlines management options for protecting wetlands during the surface mining of coal, particularly for the portion of the Eastern Interior Coal Region that is found in Kentucky, Indiana, and Illinois. It is presented in manual form for use by coal mine operators, regulatory agencies and research institutions. The previous phases of the project produced an atlas of the most heavily-mined areas of the western Kentucky coal field, which classified and identified wetlands in these areas, and discussed some specific impacts of mining on these wetlands. The need to present information that will lead to action by coal operations and regulatory agencies to protect wetland areas, is the incentive for this report. The main issues addressed in this the manual include: basic information for identifying wetlands; wetland values, and methods used for values assessment; how coal surface mining can affect wetlands; a method for addressing wetland protection needs and some prevention and mitigation actions; reclamation alternatives, including wetland restoration and the creation of wetlands as alternative ecosystems on mined areas; and general legal and regulatory information concerning wetland protection and surface mining of coal. Information was gathered through a search of current literature and by contact with state and federal agencies, some coal mining operations, and other concerned organizations. A detailed listing of places to go for more information is included as an appendix